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Abstract 

We compute the energy levels of a 2D Hydrogen atom when a constant 
magnetic field is applied. With the help of a mixed-basis variational method, 
we calculate the energy eigenvalues of the 15, 2P~ and 3D~ levels. We 
compare the computed energy spectra with those obtained via a generalization 
of the mesh point technique as well as the shifted 1/N method. We show that 
the variational solutions present a good behavior in the weak and strong 
magnetic field regimes. 
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The study of the behavior of hydrogenlike atoms in magnetic fields has been a subject 
of considerable interest in recent years. Stimulated by recent advances in nanofabrication 
technology that have made possible to create low dimensional structures containing one or 
more electrons: quantum wells, wires, and dots, a large body of papers fl|-§| (and references 
therein) has considered the problem of two-dimensional donors in semiconductor superlat- 
tices. The most commonly studied semiconductor superlattice consists of regions of GaAs 
which act as wells for the conduction electrons separated by regions of Al x Ga\- x As which act 
as barriers. Quantum wells in the presence of strong magnetic fields can form quantum dots, 
therefore the study of donor levels in Ga>i_ x Al x As would be of help in understanding the 
role played by impurity ions in low dimensional structures. The computation of the binding 
energies as well as the study of exchange and correlation effects make it necessary to have 
at our disposal reliable analytic as well as numerical methods. In this direction, different 
approaches have been considered. Besides numeric and perturbation methods |5|-[7| , we have 
to mention the two-point Pade approximation @,[|, as well as the shifted 1/N expansion 

EH- 



The Hamiltonian describing the Coulomb interaction between a conduction electron and 
a donor impurity center when a constant magnetic field B is applied perpendicular to the 
x-y plane can be written as 
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where we have chosen the vector potential A =^-G ip in the symmetric gauge. The coupling 
constant is 7 = e 2 h 3 B /(ce 3 m* 2 ) where m* is the effective mass, e is the dielectric constant of 
the host material. L z is the angular momentum operator —ihd/d<f) with eigenvalue hm.The 
units of energy are given in terms of the effective Rydberg constant TZq = m*e 4 /2h 2 e 2 and 
the effective Bohr radius a* 
The substitution 



h e/m*e , respectively. 



imtp 



u(p) 



/27T y/p 



(2) 



reduces the Schrodinger equation Hty = E^f to the following second order ordinary differ- 
ential equation 
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It is well known that exact solutions of Eq. (0) cannot be expressed in closed form in 
terms of special functions There are analytical expressions for the energy for particular 

In this paper, with the help of a mixed-basis variational method, 
and 3-D - energy levels of a 2D (x-y plane) hydrogenic 



values of 7 and m [|i~3|-|l~5[ 
we report a calculation of the IS, 2P~ 
atom in the presence of a strong magnetic field in the z direction; we compare our results 



with those obtained using the Schwartz |1| interpolation technique; as well as the shifted 
1/N method. 

In order to apply the variational method to our problem [1 
function. Since Eq. (0) reduces to the Hydrogen atom equation when 7 



we look for a trial wave 
0, we can 
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consider as a basis for 7 << 1 the Hydrogen wave functions ^h- Since < \H \ >< E„ 
we obtain an upper bound of the energy for small values of the parameter 7. The solution 
of Eq. (|3|) when 7 = is 

Mp) = 2 M • (1/2 + 2 n % W) ) W 



where -D m „ is a normalization constant, and L(a,b,x) are the Laguerre polynomials [18 
Consequently the energy spectrum in the zero-field limit takes the form, 

Eh = ~ (1/2 + n p + \m\) 2 (5) 

Conversely, for large values of 7, a good trial basis is that of the spherical oscillator. In this 
case the solution of Eq. has the form 

u 0s c{r) = C min e-^ 2 / 4 p (|m|+1/2) ^(^, \m\ , |p 2 ) (6) 
and, in the high-field limit the energy levels are 

Eqsc = 7(2n p + \m\ + m + 1) (7) 

If we attempt to apply the variational method using only the hydrogen atom basis, 
we will obtain a good agreement with the accurate results for small values of 7, but this 
approach fails for large 7 even if we consider a basis with many. Analogous situation occurs 
with the oscillator basis, for large 7, which converges very slowly for small values of 7. J7|. 

In order to solve this problem, we propose a mixed-basis approach. The idea is to use as 
trial function, for any quantum level n, a linear combination of the form 

N 

U n = ^{CiHUiH + C N+i0 UiOsc) (8) 
i 

where N > i > n; Uiu and Uio.sc are the corresponding hydrogen and oscillator wave functions 
associated with the quantum level i ; and cm are constants to be calculated. It is 
worth noticing that our mixed basis is not orthogonal under the inner product UiUjdp. 
We proceed to minimize the expectation value (u \H\ u) with the normalization condition, 
(u\u) = 1 

After performing a variation on the basis coefficients q, we reduce our problem to that 
of solving the matrix equation 

\H\ Uj) — A (ui\uj)} Cj = (9) 

after substituting the Hamiltonian ([I]) into ([|), where the lowest value of A will be the 
energy of the level. When j < 3, we can analytically compute the energy eigenvalues A 
and eigenvectors Cj. The advantage of this approach is twofold. First, we have that the 
eigenvalues satisfy the inequality A < E and therefore we have a lower bound for our energy 
levels. Second, we obtain a relatively simple expression for the normalized eigenfunctions. 
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In this paper we choose to work with a three term mixed-variational basis. In this case 
we can solve the resulting third-order algebraic equation (||) with the help of the Cardano 
method. 

In order to compute the binding energy for the ground IS state, following the above 
proposed scheme, we have two possible trial functions 

uis = cxUish + c 2 u 1S o + C3U2SO (mixl2) (10) 

and 



Wis = CxUish + c 2 u 2 sh + c 3 u 1S o (mix21) (11) 

It is worth mentioning that our three-term bases have the same angular dependence of the 
eigenfunction to be approximated. Analogously, 

u 2 p- = c\u 2 p-u + c 2 u 2 p-o + c 3 u 3P - (mixl2), u 2 p- = Ciu 2 p-H + c 2 u 3P - H + c 3 u 2 p-o (mix21) 

(12) 

and 

u-w- = cxUzd-h + c 2 u 3D - + c 3 u AD - (mixl2), u 3D - = cxUzd-h + c 2 u 4D - H + c 3 u 3D - (rnix21) 

(13) 

The numerical computations of the energy spectra associated with Eq. (13]) will be carried 



out with the help of the Schwartz method [16| This method gives highly accurate results 



given a thoughtful choice of the reference function. For Eq. (|3|) we chose as the interpolation 
function 

f(p)=J2fm ( (14) 

with u(p) = sin[7r(p//i) 1//2 ] , r m is a zero of u(p), a m is a zero of its derivative, and h is 
the step of the quadratic mesh. The use of this scheme on Eq. (|3|) leads to an algebraic 
eigenvalue problem, giving as result a non-symmetric matrix to be diagonalized in order to 
obtain the energy values. 

Here, as an illustration of the mixed-basis method, we present two different expansions 
of the IS, 2P~, and 3D~ of the 2D Hydrogen Hamiltonian (0). We plot the energy against 
7' = 7/(7 + 1) as the horizontal scale. 
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FIG. 1. Energy of the IS state as a function of 7'. The light broken line is obtained by numerical 
methods; the dashed line corresponds to the mix21 basis, (2S, 3S Hydrogen bases and IS oscillator 
wavefunction) . The solid line is obtained by using the mixl2 basis (IS, 2S oscillator bases and the 
IS Hydrogen wavefunction). The dotted line is obtained with the help of the shifted l/N method 
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FIG. 2. The figure shows the difference between the numeric result for the 15 energy spectrum 
and the energy values computed with the help of the mixl2 basis (solid line), mix21 basis (heavy 
dashed line), and the shifted 1/N method (light dashed line). 
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FIG. 3. Energy of the 2P~ state as a function of 7' The light broken line is obtained by 
numerical methods; the dashed line corresponds to the mix21 basis, (2P~ , 3P~ Hydrogen bases 
and 2P~ oscillator wavefunction) . The solid line is obtained by using the mixl2 basis (2P~ , 3P~ 
oscillator bases and the 2P~ Hydrogen wavefunction). The dotted line is obtained with the help of 
the shifted 1/N method. 
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FIG. 4. The figure shows the difference between the numeric result for the 2P~ energy spectrum 
and the energy values computed with the help of the mixl2 basis (solid line), mix21 basis (heavy 
dashed line), and the shifted 1/N method (light dashed line). 
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FIG. 5. Energy of the 3D state as a function of 7'. The solid line is obtained by numerical 
methods; the heavy dashed line corresponds to the mix21 basis (3D~ ,4D~ Hydrogen bases and 3D~ 
oscillator wavefunction) . The light dashed line is obtained by using the mixl2 basis (3D~,4:D~ 
oscillator bases and the 3D~ Hydrogen wavefunction). The dotted line is obtained with the help 
of the shifted 1/N method. 
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FIG. 6. The figure shows the difference between the numeric result for the 3D~ energy spectrum 
and the energy values computed with the help of the mixl2 basis (solid line), mix21 basis (heavy 
dashed line), and the shifted 1/N method (light dashed line). 

It is easy to see that even for a two term mixed basis a good fit is obtained in the weak 
and strong field regimes. One of the basis gives a reasonably good fitting in the intermediate 
region. A better fit is obtained with the help of the mixl2 bases. We also have that for the 
2P~ and 3D' states the mixed-basis variational approach gives very good results. Figures 
1, 3, and 5 compare the variational results with those obtained numerically and with the 
help of the shifted 1/N method. 

It is difficult to decide on the evidence of Figs. 1, 3, and 5 which technique is better 
suited to the computation of the hydrogen energy levels. Figs. 2, 4, and 6 show that the 
shifted 1/N method always gives results below the accurate energy levels. Although Figs. 
1, 3, and 5 do not show mayor changes between the variational and shifted 1/N method, 
the difference plots Figs. 2, 4, and 6 clearly show the improved performance of the mixed 
(mixl2) variational method, especially for large values of 7. 

It would be interesting to apply the mixed-basis technique for the 2D Hydrogen prob- 
lem when relativistic effects are not negligible. This will be the object of a forthcoming 
publication. 
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